Single-particle and many-body analyses of a quasi-disordered integrable system after a quench 
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In general, isolated integrable quantum systems relax to an apparent equilibrium state in which the expectation 
values of few-body observables are described by the generalized Gibbs ensemble. However, recent work has 
shown that relaxation to such a generalized statistical ensemble can be precluded by localization in the presence 
of (quasi-)disorder. Here we undertake complementary single-particle and many-body analyses of noninteract- 
ing spinless fermions and hard-core bosons within the Aubry-Andre model to gain insight into this phenomenon. 
Our investigations span both the localized and delocalized regimes of the quasi-disordered system, as well as 
the critical point separating the two. Considering first the case of spinless fermions, we study the dynamics 
of the momentum distribution function and characterize the effects of real-space and momentum-space local- 
ization on the relevant single-particle wavefunctions and correlation functions. We show that although some 
observables do not relax in the delocalized and localized regimes, the observables that do relax in these regimes 
do so in a manner consistent with a recently proposed Gaussian equilibration scenario, whereas relaxation at the 
critical point has a more exotic character. We also construct various statistical ensembles from the many-body 
eigenstates of the fermionic and bosonic Hamiltonians and study the effect of localization on their properties. 

PACS numbers: 05.30.-d, 05.70.Ln, 02.30.1k, 05.30.Jp, 05.30.Fk 



I. INTRODUCTION 

In recent years, there has been a dramatic growth in inter- 
est in the physics of non-equilibrium quantum systems, driven 
in large part by advances in experimental atomic physics, in 
particular in the area of optical lattices yj 0]. The precision 
time-dependent control and observation of quantum effects af- 
forded by these experiments, together with the high degree of 
isolation of the system from the environment, have invigo- 
rated the theoretical study of the time evolution of isolated 
many-body quantum systems. The predictions of previously 
abstract lines of theoretical inquiry into the mechanisms by 
which thermal behavior emerges from purely unitary dynam- 
ics, and the role of conservation laws and integrability in such 
dynamics, can now for the first time be directly compared with 
empirical evidence acquired in experimental laboratories. 

Several theoretical studies into the quantum origins of ther- 
malization in isolated nonintegrable systems have found that 
away from the edges of the spectrum, few-body observables 
relax to the predictions of conventional statistical ensem- 
bles O UlQIl . which has been connected Illll[l2ll to the emer- 
gence of quantum chaos 013l4l8ll . In addition, it is now well 
established that the highly constrained dynamics of integrable 
systems can in fact give rise to the relaxation of few -body ob- 
servables, and that the equilibrium values of these quantities 
are in many cases described by the generalized Gibbs ensem- 
ble (GGE) EMI ■ 

The GGE is constructed by maximizing the many-body en- 
tropy Bo, 13711 while constraining the mean values of all inte- 
grals of motion to their expectation values in the initial state 
\ipi) of the system. The density matrix in the GGE takes a 
Gaussian form similar to that of the grand-canonical ensem- 
ble, and can be written as 11911 
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where -Zqge = Tr[e ^ s ^-^j is the GGE partition function, 
the I s are the conserved integrals of motion, and the A s are the 
corresponding Lagrange multipliers, which are determined by 
the constraints Tr[p GGE I 8 ] = (*j|J g |#j). 

The significance of the fact that the GGE provides an ac- 
curate description of observables following relaxation can be 
seen by contrasting its predictions with those of the "diagonal 
ensemble" (DE) I5[]: For any initial state \^i), the time evolu- 
tion of an observable O under a time-independent (integrable 
or non-integrable) Hamiltonian H can be written as 

0(t) = J2 \C a \ 2 O aa + £ C* a C^ E ^ E ^l h O aP , (2) 
where 0{t) = (#(t)|6|#(t)>, |*(t)) 

E a Cae- US - T /'Va>. fltya) = E a \i> a ), C a = (^ a |*/>, 

and O a f} = (ip a \0\ipp). The infinite-time average of the 
observable is therefore given by 

0{t)= lim -1 f dTO{r) 
t'^oo t' J q 

= ^|C q | 2 O q « = (0)de, (3) 

a 

which defines the expectation value of O in the DE J5[]. The 
DE involves as many constraints as the dimension of the 
many -body Hilbert space (the overlaps of the initial state with 
the eigenstates of H ), which grows exponentially with system 
size. By contrast, for models that can be mapped to nonin- 
teracting Hamiltonians, the GGE involves a number of con- 
straints that is only polynomially large in the size of the sys- 
tem IU9I1 . It may, therefore, appear surprising that the predic- 
tions of the GGE for expectation values of observables can 
agree with those of the DE. 



From a many -body perspective, the success of the GGE can 
be understood as follows [28]: The eigenstates of a given in- 
tegrable Hamiltonian with similar distributions of conserved 
quantities have similar expectation values of few-body ob- 
servables (with the differences vanishing in the thermody- 
namic limit). Furthermore, the majority of the states that 
contribute to the DE have a distribution of conserved quan- 
tities similar (in a coarse-grained sense) to that of the initial 
state, and this is also the case for the states that contribute 
most strongly to the GGE. These facts imply that differences 
between the weights in the DE and the GGE are irrelevant 
and both ensembles will produce the same results for few- 
body observables in the thermodynamic limit Il28ll . This sce- 
nario can be viewed as a generalization of the eigenstate ther- 
malization hypothesis (ETH) (2H5[], and has been explored in 
Ref. IBal for describing observables after relaxation by means 
of a single representative state. 

Interestingly, it has been recently shown that in an inte- 
grable system, in the presence of localization, the GGE can 
fail to describe observables after relaxation I133ll . This effect, 
which parallels the breakdown of eigenstate thermalization in 
nonintegrable disordered lattice systems in the presence of 
localization [39], has been related to the localized behavior 
of the underlying system of noninteracting particles to which 
some integrable models can be mapped [3Qt |34j, |35|[. Here 
we gain further insights into this phenomenon by undertak- 
ing single-particle and many-body analyses of noninteracting 
spinless fermions and hard-core bosons. We study the dynam- 
ics of noninteracting fermions within the Aubry-Andre model 
previously studied in Ref. 1331 for hard-core bosons, and show 
that although the fermion momentum distribution equilibrates 
in the localized regime, it fails to equilibrate in the delocal- 
ized one. This should be contrasted with the density profiles, 
which exhibit the opposite behavior, equilibrating in the delo- 
calized regime, but not in the localized regime 13311 . We find 
that, whenever observables do equilibrate to their GGE values, 
they do so in a manner consistent with the Gaussian equilibra- 
tion scenario of Ref. J40I1 . Furthermore, we relate the failure 
of a given quantity to equilibrate in a certain regime to the 
behavior of the single-particle wavefunctions, as discussed in 
Refs. HHUl. 



For hard-core bosons, we connect the results of Ref. 113 311 . in 
which the GGE was shown to fail to describe the momentum 
distribution function after relaxation in the localized regime, 
to the single-particle results for fermions. In addition, for 
both hard-core bosons and spinless fermions, we study the 
density profiles and momentum distribution functions in the 
many-body eigenstates of the appropriate Hamiltonians. We 
focus in particular on the behavior of these quantities in the 
many-body eigenstates that contribute to the DE, to the mi- 
crocanonical ensemble (ME), and to the microcanonical ver- 
sion of the GGE; i.e., the generalized microcanonical ensem- 
ble (GME) |28|[ . We find indications that single-particle real 
space localization in the localized phase and momentum space 
localization in the delocalized phase lead to a distinctive be- 
havior of the many-body eigenstate expectation values of the 
density and the momentum distribution of the fermions, re- 
spectively. No similar effect is detected in the many-body 



eigenstate expectation values of the momentum distribution 
of the hard-core bosons. 

This paper is organized as follows. In Sec. HU we introduce 
the models, quench protocols, and observables studied in later 
sections. We also review the statistical ensembles utilized to 
describe observables after relaxation. The time evolution of 
spinless fermions following a quench is studied in Sec. [HI] 
Specifically, we examine the relaxation dynamics and time 
fluctuations of one-body observables, as well as properties of 
the single-particle eigenstates that help us understand the ob- 
served out-of-equilibrium behavior. Section [TV] is devoted to 
the study of one-particle observables in the many -body eigen- 
states of the bosonic and fermionic Hamiltonians. In SeclVl 
we summarize our results and present our conclusions. 



II. MODELS AND QUENCHES 

We consider two models on a one-dimensional lattice with 
open boundary conditions: noninteracting spinless fermions 
(SFs) and hard-core bosons (HCBs). The HCB model can 
be mapped onto the model of SFs, by mapping it first onto a 
spin- 1/2 chain via the Holstein-Primakoff transformation Il4jll 
and then onto SFs via the Jordan- Wigner transformation I142ll . 
In both cases, we study the effects of an additional periodic 
potential, with a period incommensurate with that of the un- 
derlying lattice of the tight-binding model, which results in 
the well known Aubry-Andre model B43I1 . The Hamiltonians 
for SFs and HCBs are given by 

L-X L 

H f = -tJ2 (f}fj+i + H - c -) + A E cos(2tt^>J, (4) 
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respectively, where L is the length of the lattice, fj and /T 

(L and by) are fermionic (bosonic) annihilation and creation 

operators on site j, and h- = nfj (n b j = b\bj) are SF (HCB) 
site occupation number operators. The prohibition of multiple 
occupancy of a single site for HCBs is enforced by the hard- 
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0. We denote the hopping param- 



eter by t, and the strength of the incommensurate potential by 
A. To enforce the incommensurability of the lattice potential, 
we use an irrational value for <;. We select the inverse golden 
mean, ? = (y/E— l)/2, which is considered to be the most ir- 
rational number B44II . For each given system size, we take the 
total number of SFs (HCBs) N f (N b ) to be N f = N b = N. 
In what follows we set t to unity; i.e., we take t as our unit 
of energy, and we also set h = 1 and the Boltzmann constant 

The fermionic Hamiltonian <j4j is quadratic and there- 
fore trivially solvable: All many-body eigenstates can be 
constructed as Slater determinants of the single-particle en- 
ergy eigenstates in the incommensurate periodic potentials. 



Although the bosonic Hamiltonian (0 is also superficially 
quadratic, the hard-core constraints on the bosonic cre- 
ation and annihilation operators encode interactions between 
the bosons, precluding a direct diagonalization in terms of 
single-particle states. Nevertheless, it can be solved via the 
combined Holstein-Primakoff and Jordan- Wigner transforma- 
tions, which implies that SF and HCB systems with Hamilto- 
nians Eqs. (0]i and 01, respectively, share the same (many- 
body) energy spectrum and consequently have identical ther- 
modynamic properties. Moreover, the two models have iden- 
tical site occupations h- = fih = hj. 

The properties of HCBs in the Aubry-Andre model have 
previously been investigated, at both zero 11451 14611 and finite 
temperature 114711 . The single-particle Aubry-Andre model 
exhibits a localization-delocalization transition at the criti- 
cal potential strength A c = 2 IJ430 : All single-particle states 
are extended when A < A c and exponentially localized when 
A > A c . At the critical point, the energy spectrum exhibits 
the fractal structure of a Hofstadter butterfly lEJJll . As HCBs 
can be mapped onto non-interacting fermions, they of course 
inherit this phase transition when subjected to the incommen- 
surate Aubry-Andre potential. In the localized regime, cor- 
relations in the ground state of the bosonic system decay ex- 
ponentially with spatial separation, and the system is said to 
form a Bose glass y|]. 

Here we study the dynamics and behavior of observables 
following sudden quenches of the incommensurate lattice 
strength A. We select initial states \^> i) to be the ground 
states of initial Hamiltonians Hi with parameters Xj, and 
the dynamics are followed under final Hamiltonians Hp with 
parameters Xf, corresponding to an instantaneous change 
(quench) of the lattice strength from Xj to Xp. We focus on 
one-body observables that are accessible in optical lattice ex- 
periments: the density profiles fij, and the momentum distri- 
bution functions 



Diagonal ensemble. The density matrix of the DE is defined 
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We note that although the site occupations are equal for SFs 
and HCBs, the off-diagonal spatial correlations, and therefore 
the momentum distributions, of the two systems are distinct. 
To calculate (wjfc) in a pure state, we follow the approach 
of Refs. 1149145111 . while for calculations in GGE we use the 
methodology of Ref. 15211 . 



Ensembles of interest 

We characterize the behavior of our quasi-disordered sys- 
tem after relaxation following a quench by comparing it to 
three different statistical ensembles: DE, ME, and GME, each 
of which we briefly describe here. 
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(8) 
i.e., it is diagonal in the (many-body) energy representation, 
and the energy eigenstates are weighted according to their 
overlaps with the initial state. The expectation value of an 
observable in this ensemble is given by Eq. 01. In various 
computational studies 1161171. 12811 . observables after relaxation 
have been shown to approach the DE predictions [Eq. (0] as 
the system size increases. 

Microcanonical ensemble. The density matrix of the ME 
can be written as 
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i.e., all eigenstates in the energy window [E — <5ME,-Eini + 
Sme] (°f which there are AFeJme) are given equal weight. We 
have checked that expectation values of observables within 
our microcanonical calculations are robust against the exact 
value of Sme- To ensure this, we select Sme to be much smaller 
than the full spectrum width, but sufficiently large to contain 
many eigenstates (i.e., larger than the average level spacing 
at the given E). In general, Sme = 0.05 for most results 
reported in this work. 

Generalized Gibbs ensemble. The density matrix for this 
ensemble, which is of a similar Gaussian form to that of the 
grand-canonical ensemble, was already introduced in Eq. (0. 
We note that a recipe for constructing the appropriate con- 
served quantities, allowing for the extension of the GGE de- 
scription to more general systems, has recently been pro- 
posed 1530 . However, we make here the 'natural' choice for 
the conserved quantities, taking them to be the occupations of 
the single-particle eigenstates of the noninteracting SFs [19]. 

Generalized microcanonical ensemble. The GME is the 
microcanonical version of the GGE 02811 . The only energy 
eigenstates that contribute to this ensemble are those that have 
distributions of the conserved quantities that are similar (in a 
coarse-grained sense) to that of the initial state. These eigen- 
states are all assigned the same weight, as in the usual ME. 
The density matrix of the GME therefore has the form 
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where S a measures the distance of the eigenstate \t/j a ) from 
the target distribution of conserved quantities determined by 
the initial state, and.A/{j s } i( s GME is the number of energy eigen- 
states within the GME window S a < Sgme- 

In order to construct the GME, one needs to compare the 
distribution of conserved quantities in each of the eigenstates 
of Hp with that of the initial state. Since the conserved quan- 
tities are fermion occupation numbers of single-particle en- 
ergy eigenstates, their expectation values in the many-body 



eigenstates I s<a = (ip ] a \ I s \i/j ' a ) are either or 1. By contrast, 
the occupations of the conserved quantities in the initial state 
(same as in the DE) can assume any real value between these 
two values; i.e., < (/ s )de = J2 a |Cq| 2 ^,q < 1. To com- 
pare those distributions, in a coarse grained way, we proceed 
as follows (H: 

(i) We sort the conserved quantities in such a way that 
(I S ) DE decreases monotonically as s increases. 

(ii) After sorting, we generate a discrete target distribution 
of conserved quantities from (J s )de- For this, we interpo- 
late (is) de t° find a continuous function I(s) (where s can 
now be any real number in the interval [0.5, L + 0.5]), and 

then compute all values s ; * satisfying L l , I(s) ds = 0.5 and 

/ » I(s) ds — 1 for I > 1. The set {.s^}, together with the 

set of corresponding weights {/(s*)}, defines the target dis- 
tribution. 

(iii) We introduce a measure 8 a to quantify how close the 
distribution of each many -body eigenstate \ip a ) is to the target 
distribution. Those states with S a < 8gme are included in the 
GME. The choice of S a is not unique. As in Ref. 12811 . we 
consider 



N 



1 



N 

TV ^ 
i=i 



I(st)(\* 



(11) 



where s t , with I = 1, . . . , N, is the number of the sorted 
single-particle states [see (i)] occupied in eigenstate \ip a }- The 
value chosen for (Sgme is the one that leads to the minimum 
value of the normalized absolute difference 
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(12) 



between the expectation values of the conserved quantities in 
the GME and DE. Similarly to the microcanonical calcula- 
tions, we have checked that our results are robust against small 
changes of Some- 



III. SINGLE-PARTICLE ANALYSIS 



In Ref. 13311 it was shown that, following a quench of HCBs 
to the localized regime of the Aubry-Andre model, one-body 
observables that depend on nonlocal correlations, such as m b k , 
do relax to time-independent values (with fluctuations van- 
ishing in the thermodynamic limit), but these values are not 
consistent with the predictions of the GGE. This should be 
contrasted with the on-site density, whose time average agrees 
with the GGE results for all regimes. It was also found that the 
dynamics of n,j and m b k are qualitatively different in the de- 
localized and localized regimes: The momentum distribution 
m\ approaches a time-independent value with increasing sys- 
tem size regardless of whether the system is in the delocalized 
or localized regime, while rij only exhibits such relaxation in 
the delocalized regime. All results for the density also apply 
to SFs, to which HCBs can be mapped. 



Here, we begin by studying the relaxation dynamics of the 
momentum distribution m k of SFs, which is in general com- 
pletely unrelated to the momentum distribution of the corre- 
sponding system of HCBs. In fact, the GGE describes, by con- 
struction, the infinite-time averages of all one-body fermionic 
observables regardless of whether the single-particle states 
are localized or delocalized and independently of the sys- 
tem size. This can be straightforwardly proven by projecting 
/5(t) = |^> (t)) (\& (t)| onto the single-particle sector. Consid- 
ering that all eigenstates of the many-body Hamiltonian are 
(antisymmetrized) direct products of the single particle states 
\s) in which Hf is diagonal (Hf\s) = e s \s)), the time evolu- 
tion of the one-particle density matrix can be cast in the form 
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In the absence of degeneracies in the single-particle spec- 
trum (which is the case in the Aubry-Andre model considered 
here), the infinite-time average of /5 sp (r) can be written as 



p sp (r) = lim — / dr, 



*pM = y>~i*x*i, d4) 






which is, by construction, the single-particle density ma- 
trix predicted by the GGE, as c ss = J2 a \C a \ 2 Is,a = 
Tr[/5cGE is]. We emphasize, however, that this does not nec- 
essarily imply that all such observables exhibit relaxation to 
their GGE values. In particular, the results of Ref. IB3II in- 
dicate that the on-site densities in the localized phase exhibit 
finite fluctuations about their GGE expectation values even in 
the thermodynamic limit. 

In order to quantify how closely ?ti{(t) = 

($(r)|m^|^(r)) approaches the corresponding GGE 
prediction, we compute the normalized difference 
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between the instantaneous momentum distribution at time r, 
and the GGE prediction for this quantity. Relaxation of m k 
to the GGE prediction is observed if 8mf (r) vanishes at long 
times, in the thermodynamic limit. In finite-sized numerical 
calculations, however, 8m' {t) will always be finite, because 
of finite-size effects. The signature of relaxation to the GGE 
in our calculations is therefore that 5m' (t) fluctuates about a 
finite average value at long times, and that this average value 
scales towards zero with increasing system size. 

In Fig. [TJ we show results for 8m' (r) in quenches with 
Xj = 0; i.e., a delocalized initial state (left panels), and 
Xj = 8; i.e., a localized initial state (right panels). After the 
quench, Xp = 1 [delocalized regime, Figs. [TJa) andQJe)], 
Xp = 2 [critical point, Figs.[TJb) and[TJf)], and Xf = 3 and 
4 [localized regime, Figs. []I C X EgX and [TJd), [TJh), respec- 
tively]. The results presented correspond to three different 
system sizes (L = 10, 100, and 1000, curves from top to bot- 
tom in each panel, and N — L/2), and to the same quenches 
and system sizes studied for HCB systems in Ref. 113311 - 
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FIG. 2: (Color online) Finite-size scaling of <5m/(oo) (see text) for 
the quenches studied in Fig.Q] Thin continuous lines depict power- 



law fits in which 5m* (oo) oc L 



for quenches with Xf 



X c and Sm f (oc) oc £-°-5°°±°-°°s in quenches with \ F > \ c . The 
power laws were fitted to results with L in the interval [200, 1000]. 



FIG. 1 : (Color online) Relaxation dynamics of m k in quenches A = 
Xi — ► Xf (as indicated in the panels) for systems with 10, 100, 
and 1000 lattice sites (curves from top to bottom in each panel), and 

N = L/2. 



The results obtained for a given final value of the incom- 
mensurate potential strength are qualitatively similar, inde- 
pendently of whether the initial state is delocalized or local- 
ized. In quenches to the localized regime [Figs. [He), (d), (g), 
(h)], we observe that 5m' (r) decays to a finite value, about 
which it undergoes fluctuations, and that this value decreases 
with increasing system size. In quenches to the delocalized 
regime [Figs. 02a) and[TJe)], 8m' {t) similarly undergoes de- 
cay to a finite value about which it fluctuates. However, in 
the delocalized case, the value to which 5m* (r) decays does 
not appear to exhibit such a pronounced reduction as the sys- 
tem size is increased, suggesting that it may not tend towards 
zero as L —$■ oo. Following quenches to the critical point 
[Figs-Qlt)) andQJf)], 5m* (t) exhibits behavior similar to that 
observed in the localized regime, decaying to exhibit fluctu- 
ations about a constant value which decreases with increas- 
ing system size. However, in this critical regime the decay of 
5mf (t) is much slower and, e.g., fluctuation about a constant 
value in the case L = 1000 is only obtained for r > 10 5 . 

To gain a quantitative understanding of the dependence 
of the long-time behavior of Eq. ( fl~5T > on the system size, 
we consider the average of 5m* (t) over the time interval 
r G [10 5 , 10 6 ], which we take as representative of the aver- 



age value after any initial transient dynamics have subsided 
— such that 5m* (r) is simply fluctuating about a constant 
value — and away from any revival, and denote this quantity 
by 5m* (oo). In Fig. [2 we plot 5m* (oo) against L for all the 
quenches shown in Fig. [T] The scalings make apparent that 
<5m/(oo) converges to a finite value as the system size is in- 
creased in quenches to the delocalized phase; i.e., although 
the infinite time average of each m^ trivially agrees with its 
expectation value in the GGE, its instantaneous value does 
not relax to the GGE in the thermodynamic limit. By con- 
trast, in quenches to the localized regime 5m* (oo) decreases 
with increasing system size as a power-law, which we find to 
be close to L~° ,50 . In the quenches to the critical regime, 
5m* (oo) is much slower to reach a clear power law scal- 
ing, but for large system sizes its behavior is consistent with 
5m* '(oo) ex L-°- 25 . 

The behavior of 5m* (oo) in the delocalized and localized 
regimes is exactly opposite to that observed for the density 
in Ref. IB3J p4ll . There it was found that the the normal- 
ized difference of the density 5n(oo) oc L~ - 5 in quenches 
to the delocalized phase, whereas it converges to a finite value 
in quenches to the localized phase. Intuitively, the extended 
states in the delocalized regime of the quasi-disordered lattice 
can be thought of as states that are localized in momentum 
space. In fact, the correspondence between position (momen- 
tum) space results in the delocalized phase and momentum 
(position) space quantities in the localized phase can, in the 
case of free fermions, be seen to be an exact consequence of 
the well-known self-duality of the Aubry-Andre model 114411 . 
We note also that the behavior of 5m* (oo) at the critical point 
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FIG. 3: (Color online) Histograms of the time fluctuations of rij (left 
panels) and m k (right panels) in quenches with Xi = and \f = 1 
(a),(d), Xf — 2 (b),(e), and Xf — 4 (c),(f). In all panels, we report 
results for systems with L = 200 and L — 1000, and Gaussian 
fits to the data for L = 1000. From the fits, the means were found 
to be zero (±0.003) in all cases and the variances were found to be 
a = (6.770 ± 0.004) x 10" 3 (a), (3.206 ± 0.009) x 10~ 2 (b), 
(1.008 ± 0.007) x 10 _1 (c), (2.25 ± 0.02) x 10~ 2 (d), (2.322 ± 
0.005) x 10 _1 (e), and (1.336 ± 0.001) x 10" 2 (f). 



FIG. 4: (Color online) Histograms of the time fluctuations of rij (left 
panels) and m k (right panels) in quenches with Ai = 8 and Xf — 1 
(a),(d), Xf = 2 (b),(e), and Xf — 4 (c),(f). In all panels, we report 
results for systems with L — 200 and L = 1000, and Gaussian 
fits to the data for L = 1000. From the fits, the means were found 
to be zero (±0.003) in all cases and the variances were found to be 
a = (1.2373 ± 0.0004) x 10~ 2 (a), (2.570 ± 0.008) x 10~ 2 (b), 
(2.27 ± 0.01) x 10" 2 (c), (4.318 ± 0.005) x 10" 2 (d), (1.860 ± 
0.004) x 10 -2 (e), and (3.360 ± 0.003) x 10" 3 (f). 



is similar to that found for <5n(oo) in Ref. 13311 . 

Aside from the peculiar behavior observed at the critical 
point, which is understandable given the very special charac- 
ter of the single-particle problem in this regime, it is remark- 
able that whenever relaxation of rij or mf, takes place (in the 
delocalized or localized regime) the time fluctuations are pro- 
portional to L~ - 5 . This is consistent with the Gaussian equi- 
libration scenario of Ref. f|40j], in which the square root of the 
normalized time variance of one-body observables in nonin- 
teracting fermion systems was argued to scale as 1/vZ 115511 . 

In Figs. [3] and |U we present histograms of the distributions 



of differences Srij (r) 



,(r) 



>GGE between the in- 



stantaneous values of the site occupations and their mean val- 
ues in the GGE (left panels), and the analogous quantities 



,/ 



(r) 



,/■ 



gge calculated for momentum- 



6m{(T) 

mode occupations (right panels). These histograms represent 
the full distribution of fluctuations of the occupations of all 
lattice sites j, and all momentum modes k, in quenches with 
A/ = (Fig.O and A/ = 8 (Fig.[4]i. Once again, the results 
for a given value of Xf can be seen to be qualitatively similar 
independently of the initial state. The histograms of the time 
fluctuations of rij in quenches to the delocalized regime and 
of m k in quenches to the localized regime have a Gaussian 
shape with a width that decreases with increasing system size. 
By contrast, the histograms of the time fluctuations of rij in 
quenches to the localized regime and of m^ in quenches to 
the delocalized regime are in general non-Gaussian (as can be 
seen by comparing them to the indicated best Gaussian fits 



to the distributions) and the widths of the distributions are 
not seen to decrease with increasing system size. We note 
that the distributions of time fluctuations of individual lattice- 
site (momentum-mode) occupations in the localized (delocal- 
ized) phase, which we have not shown, are quite strongly non- 
Gaussian, exhibiting, e.g., bimodal structures. The behavior 
of the time fluctuations at the critical point is intermediate be- 
tween what is seen in the localized and delocalized phases: 
The fluctuations of rij and m£ are close to Gaussian, with a 
width that decreases less dramatically with increasing system 
size, consistent with the L~ - 25 behavior found for Sin' (oo) 
(Fig. |2). This exotic behavior at the critical point warrants 
a more specific investigation that is beyond the scope of this 
article. 

An understanding of why Gaussian equilibration fails to oc- 
cur for both fermionic observables at the critical point, for rrr k 
in the delocalized phase, and for rij in the localized one, can 
be gained through an analysis of the properties of the single- 



particle eigenstates of Ht 
mentum space OCX [34 l35l 



\s) = 7]|0)) in both real and mo- 
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^fescjj (cl creates a SF at momentum fc), it follows that 
: (i) if \s) is delocalized in real (momentum) space then 
| 2 ~ 1/L (|Wfe«| 2 ~ 1/L) and a 2 ,, (a^ ) must decrease 
as 1/L or faster (because ^ s , \p I ss , | 2 = Nf) with increas- 
ing system size, and (ii) if | s) is localized in real (momentum) 
space then the the corresponding variance er 2 . (o- 2 ^ ) will be 
dominated by rare large values of \p ss , | 2 and will remain finite 
in the thermodynamic limit. 

In the insets to Fig. |5j we plot the average of the square 
root of the variance over all lattice sites <r* ve = 1/L J\. a nj 
[insets in Fig. |3a)] and over all momentum states cr™ = 
l/LJ2k a m k [insets in Fig.|5Jb)] against the system size L, 
for systems with between 10 and 1000 lattice sites and sub- 
ject to the same quenches studied previously (A/ = in the 
top insets and Xj = 8 in the bottom ones). The results for 
a™ e and a^ can be seen to be qualitatively similar to those 



in Ref. 13311 for 8n{oo) and in Fig. [2] for 6m* (oo). In order 
to relate the behavior of the variances to the properties of the 
single-particle eigenstates \s) in real and momentum space, 
we compute the average inverse participation ratios (IPRs) 



IPR » = iEE 

* 3 



IPR 



\Vks 



(18) 



s fc 



The results for the IPRs are reported in the main panels in 
Fig. [5] They show that, as expected, the eigenstates of the 
Hamiltonian are delocalized in real space (IPR„ ~ 1/L) and 
localized in momentum space (IPR m ~ LP) for Xf < 2, 
and localized in real space (IPR„ ~ LP) and delocalized in 
momentum space (IPR m ~ 1/L) for Xp > 2, which ex- 
plains the behavior of the time fluctuations of the density and 
fermionic momentum distributions in the delocalized and lo- 
calized regimes. For Xp = 2, we find that IPR n ~ IPR rn ~ 
1/y/L, illustrating the exotic structure of the eigenstates of 
the critical Hamiltonian in both real and momentum space. 
From this scaling we can infer that \uj s \ 2 ~ \vk s \ 2 ~ 1/L - 75 
when Xp = 2, implying that ct 2 . and er 2 nfc decay like l/s/L 
or faster at the critical point, which is indeed consistent with 
the results of Fig. [2] We note also that in general, the aver- 
age square-root variances and IPRs in the critical regime reach 
their exotic scaling limits at quite small values of L, compared 
to the behavior of 6m* (oo) at the critical point (Fig. |2J. The 
unambiguous scaling at criticality seen in Fig. [5] thus lends 
strong support to our identification of 6m* (oo) as scaling like 
~ L-° 25 . 

For fermionic observables, we should stress the fact that, in 
the localized regime, localization in real space precludes re- 
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FIG. 5: (Color online) Inverse participation ratios in real (a) and mo- 
mentum space (b) in the delocalized and localized phases as well as at 
the critical point. Thin continuous lines are power-law fits in which 
(a) IPR„ oc L- lh99±0 - 02 for \ F = 1 and IPR n oc L-°- 52±om 
for \ F = 2, and (b) IPR m oc £-°-6B±o.oi for Xp = 2 and 
IPR„ oc ^-i-OOBiO.ooi for Af , = 4 The insets in (a) depict jm 

for quenches with A/ = (top inset) and Xi = 8 (bottom inset), 
while the insets in (b) depict <j^ e for quenches with A/ = (top in- 
set) and Xj — 8 (bottom inset). Thin continuous lines are power-law 
fits in which, in both insets in (a), aZ oc £-°-«>±o.oi for Xp = 1 
and a'T oc L^ a25±0 ' 01 for A F = 2, and, in both insets in (b), 



-0.500±0.005 



aZ oc L-"-^±u. UJ for Xf = 2 and o'Z oc L-^-^uztu.uuo for 
Af = 4. The power laws were fitted to results with L in the interval 
[160, 1000]. 



laxation of the density profiles in the same way that, in the de- 
localized regime, localization in momentum space precludes 
relaxation of the momentum distribution. This symmetry be- 
tween localization in real and momentum space is broken in 
the case of HCBs, because the mapping to the underlying 
model of SFs only preserves correlations that are diagonal in 
real space (i.e., properties related to the density). Thus for 
HCBs, although relaxation of the density is precluded in the 
localized regime, relaxation of the momentum distribution can 
occur in the delocalized regime, as was found in Ref. 0311 . 

As discussed in Refs. 11301,1341, 13511 . if the variances of one- 
particle correlations (here we have focused only on the den- 
sity) do not vanish with increasing system size — which is 



only possible if off-diagonal elements of \p I ss , | 2 contribute tc 
the fluctuations in the thermodynamic limit — then Wick'; 
theorem can break down for time averages of higher order cor- 
relations. We recall that nonlocal one-particle correlations o: 
HCBs [e.g., Eq. ©J correspond, via the Jordan-Wigner trans 
formation, to higher order correlations of SFs [13511 . Thus the 
breakdown of the GGE for describing m h k after relaxation ir 



the localized phase of HCBs 0311 can be understood as a direc 
consequence of the fact that time fluctuations of local one 
particle correlations of the underlying free-fermion model re 
main finite as one increases the system size in the localizec 
phase. 



IV. MANY-BODY ANALYSIS 

In this section, we study the density profiles and momenturr 
distribution functions in the many-body eigenstates of the SI 
and HCB Hamiltonians. Our goal is to understand how local- 
ization, in real and momentum space, affects the many-body 
eigenstate expectation values of these observables. Since this 
study requires the construction of the full set of energy eigen- 
states of the many-body system, which grows exponentially 
with increasing system size, our analysis will be restricted to 
lattice lengths that are much smaller than those studied in the 
single-particle analysis of the previous section. The small- 
est systems considered in this section have 20 sites and the 
largest ones have L = 50, with N = L/5, as opposed to the 
filling N = L/2 (which yields the maximal Hilbert space di- 
mension for a given system size) considered in the previous 
section. The largest systems we consider have a Hilbert space 
of dimension O(10 10 ). In order to compare systems that have 
equivalent excitation energies per particle after the quench, we 
will focus on three quenches that, while having the same final 
Hamiltonians as in the previous section, have for their ini- 
tial states the respective ground states of Hamiltonians with 
A/ = -1.5 for X F = 1.0, A/ = -0.5 for X F = 2.0 and 
\j = 1.0 for Xp = 4.0. These three quenches lead to time- 
evolving states whose energies are similar to those of systems 
in thermal equilibrium with T ~ 1.7. 

In Fig. [6] we show the distribution of conserved quantities 
in the quenches mentioned above for systems with L = 50. 
We compare results for those distributions in the initial state 
(same as the DE and the GGE), in the GME (constructed as 
described in Sec. [H), and in the ME. The contrast between 
the distribution of conserved quantities in the initial state and 
in the ME is apparent, whereas the distribution in the GME 
closely agrees with that in the initial state. Figures |6jd)- 
|6ji) depict how the normalized absolute differences between 
the distribution of conserved quantities in the GME and ME 
and that in the initial state [Eq. ([T2l and the analogous ex- 
pression obtained by replacing (I s )gme with (J s )me, respec- 
tively] scale with increasing system size. For the systems 
analyzed here, these differences are always smaller for the 
GME than for the ME. More importantly, for all quenches, 
they are seen to decrease with increasing system size for the 
GME [Figs.|6]:d),Hf), andith)]. The differences between the 
ME and the DE exhibit clear saturation behavior in the delo- 
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FIG. 6: (Color online.) [(a),(b),(c)] Conserved quantities (I s ) in 
the DE, GME, and ME for L = 50, N = 10, sorted in order of 
decreasing occupation in the DE. [(d)-(i)] Scaling of the absolute 
differences D between distributions of conserved quantities with in- 
creasing system size. [(d),(f),(h)] Absolute differences D between 
the conserved quantities in the GME and the DE [see Eq. <U2H. 
[(e),(g),(i)] D between the conserved quantities in the ME and the 
DE [same as Eq. i ll 2b but with "GME"— > "ME"]. The quenches are 
indicated as Xi — > Xf- For quenches A/ = —1.5 — > Xf — 1.0 and 
Xi — —0.5 — > Xf — 2.0 we found that <5gme = 0.85 yields the 
minimum value of D, and for quenches Aj = 1.0 — > Xf = 4.0 we 
found Some = 0.95. For the ME, <5me = 0.05 in all quenches. 



calized and critical regimes [Figs. |6je), |6lg)] ■ In the localized 
regime [Figs. |6ji)], the difference between the ME and DE is 
both smaller than, and does not exhibit saturation quite as ob- 
viously as, that in the other two regimes. We note, however, 
that it was found previously in Ref. [28] that the discrepancy 
between distributions of conserved quantities in the ME and 
DE can be strongly dependent on the initial state (and in par- 
ticular on its energy). 

Although disagreement between the distributions of con- 
served quantities in the DE and ME indicates the failure of 
the ME to describe the state of the system after relaxation, 
the degree of agreement between the corresponding distribu- 
tions in the DE and GME yields only incomplete information 
on the accuracy of the GGE as a characterization of the sys- 
tem at long times. In particular, the distributions here apply 
equally to SFs and HCBs, whereas the results of Ref. IB3I1 and 
Sec. Un] indicate that the presence or absence of relaxation, 
and the agreement between time averages and the GGE pre- 
dictions, can differ between SFs and HCBs with equal quench 
parameters. To further characterize the relationship between 
(generalized) thermalization and the structure of the many- 
body eigenstates of the system, we now turn our attention to 
the behavior of the density and momentum distributions in the 
eigenstates of the SF and HCB systems. 

In order to quantify the differences between the predictions 
of the GME and ME for each observable and those of the DE, 



we compute the normalized difference, 
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where O is either the site occupation n (for which t = j) 
or the momentum occupation nv , m b (for which £ = k), 

O s taM = E Q < M (^a|^|V«) S <0,) stat , With < tat the 

weight of each many-body eigenstate in the relevant ensem- 
ble, and 'stat' stands for one of the three ensembles: the DE, 
the GME, or the ME. Hence, the agreement between a statisti- 
cal ensemble and the DE in the thermodynamic limit becomes 
apparent if AO stat ~ DE vanishes with increasing system size. 

We quantify the behavior of the eigenstate expectation val- 
ues of the observables by calculating the average variance 
within each ensemble, which we define by 
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tat,£ = l^a^aWa^eWa)" P (O'Dstat- A finite 

) <Tq E in the thermodynamic limit implies that 



(generalized) eigenstate thermalization does not occur. How- 
ever, the vanishing of (<7q ME ) (Jq E is not a sufficient condi- 
tion for (generalized) eigenstate thermalization to occur: The 
(GME) ME may still fail to describe observables after relax- 
ation if the differences between the expectation values of ob- 
servables in distinct eigenstates contributing to the ensemble 
do not vanish in the thermodynamic limit. In such a scenario, 
the results of the DE for the expectation values of observables 
may be dominated by "rare" states, and thus be distinct from 
those of the (GME) ME [8J]. We calculated the individual dif- 
ferences between observables in each eigenstate and in the 
ensemble average and attempted to quantify the scale of the 
maximal differences with increasing system size. However, 
the results were found to be dominated by finite size effects 
and we could not extract any consistent information from our 
investigations. We therefore only report results for cr^" in this 
section. 

We further note that, in a previous study of HCBs systems 
whose properties are qualitatively similar to those of the sys- 
tems studied here for A < A c 1281 . indications were found 
that c^ 18 (for m b k ) vanishes with increasing system size while 
(rh b k )GME = ("Jfe)DE- Since the weights w GME and iu° E were 
seen to be different 02811 . with an exponentially smaller num- 
ber of states usually contributing to the DE when compared to 
the GGE JH], the results of Ref. JUt] hint that a generalized 
eigenstate thermalization may be at play for m k . Namely, that 
all eigenstates with similar distributions of conserved quanti- 
ties have similar expectation values of the HCB momentum 
distribution function (with differences that vanish in the ther- 
modynamic limit). 

We now present results that indicate that this mechanism 
is indeed the explanation why the GGE provides the correct 
description of HCB observables after relaxation in quenches 
to the delocalized phase. In Fig. [7] we show density plots of 
the coarse-grained weights with which eigenstates contribute 
to the DE [Fig. Ha)] and to the GME [Fig. 13b)] for delocal- 
ized HCB systems with L = 45 (main panels) and L = 35 
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FIG. 7: (Color online.) Density plots of the coarse-grained weights 
of eigenstates in (a) the DE and (b) the GME as a function of the 
expectation value of m(k = 0) and the eigenstate energy per site 
e a = E a /L, for quenches to the delocalized regime. Results are 
shown for systems with L — 45 (main panel) and L = 35 (insets). 
The width of the windows used in the coarse graining is 0.005 in k 
and 0.002 in e. 



(insets). We see that in both ensembles, as the system size in- 
creases, weight becomes increasingly concentrated in eigen- 
states with — 1.0 < e a = E a /L < —0.5 (though more clearly 
in the GME than in the DE). Moreover, the expectation val- 
ues of rh(k = 0) in these most highly weighted eigenstates 
are narrowly distributed compared to the full range of expec - 
tation values of all eigenstates within the same energy range. 
Furthermore, it is remarkable that the expectation values of 
rh(k = 0) in the dominant states of the two ensembles are 
similar to each other, and that this agreement is seen to im- 
prove with increasing system size. 

In Fig. [8] we present the corresponding results for quenches 
to the localized phase. We observe qualitative differences 
between the distributions of weights in the DE and in the 
GME: the former are spread relatively smoothly over energies 
e a < 2 [Fig.[HJa)], whereas the latter tend to concentrate in 
several distinct energy bands [Fig.[5Jb)]. Moreover, in the DE 
the weights tend to increase as e a decreases, corresponding to 
larger values of the expectation of m(k = 0), whereas in the 
GME the weight distribution is more strongly concentrated 
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FIG. 8: (Color online.) Density plots of the coarse-grained weights 
of eigenstates in (a) the DE and (b) the GME as a function of m(k — 
0) and the eigenstate energy per site e a = E a /L, for quenches to 
the localized regime. Results are shown for systems with L = 45 
(main panel) and L — 35 (insets). The width of the windows used in 
the coarse graining is 0.001 in k and 0.004 in e. 



in eigenstates with higher values of e a , and smaller values 
of the zero-momentum expectation values. This suggests that 
the dominant eigenstates in each of the two ensembles do not 
yield similar momentum distribution functions, implying the 
failure of generalized eigenstate thermalization in the local- 
ized phase, similar to the failure of eigenstate thermalization 
in nonintegrable systems in the presence of localization II39I1 . 
In what follows, we explore in more detail the properties of 
the eigenstate expectation values of the various observables of 
interest for SFs and HCBs, and how they change with increas- 
ing system size. 



Quenches to the delocallzed regime 

We start by studying the behavior of the density profiles in 
the many-body eigenstates of the SF and HCB Hamiltonians. 
The density profiles of SFs and HCBs are identical, and, by 
construction, the predictions of the GGE for expectation value 
of this observable are the same as those of the DE. 

In Fig. |9|a), we show the density profiles after a quench 
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FIG. 9: (Color online.) (a) Density profiles (L = 50), (b) An* DB , 
and (c) o-jjf 1 in quenches to the delocalized regime, \i = —1.5 — ¥ 
Xf = 1.0. In (b), the main panel depicts An GMME , while the inset 
depicts An MME . The legends for the ensembles considered are the 
same in all panels and are reported in (a). The values of Some and 
<5me are the same for all L's and are equal to the ones in Fig. [6] 



to the delocalized regime as predicted by the DE, the GME, 
and the ME. The agreement between DE and GME is excel- 
lent, leading to the small values of An GMEDE in Fig. |Hb). 
The latter quantity is seen to decrease with increasing system 
size indicating that the DE and GME predictions will agree 
in the thermodynamic limit. On the other hand, in Fig. |9}a), 
large differences can be seen between the outcomes of DE 
and ME calculations for the density profiles, which leads to 
large values of An MEDE as depicted in the inset in Fig.|9jb). 
The results in the inset suggest that An ME " DE saturates to a fi- 
nite value with increasing L, which would be consistent with 
the results of Ref. 13311 . where the grand-canonical ensemble 
was shown to fail to predict the outcome of the relaxation dy- 
namics with increasing system size for this observable (for 
systems up to 20 times larger than those considered here). 

Figure[9]c) shows the scaling of the average variance of the 
site occupations in the DE, the GME, and the ME with L. In 
all ensembles, this variance is of a similar small magnitude, 
and decreases with increasing system size. However, our re- 
sults are not conclusive as to whether the variance (in any of 
the ensembles) vanishes or saturates to a finite value in the 
thermodynamic limit. Inasmuch as we understand thermal- 
ization in an isolated system to result from eigenstate thermal- 
ization, the fact that this observable does not thermalize 113311 
implies that eigenstate thermalization does not occur in this 
system. Given the fact that the predictions of the GGE and 
GME for this observable agree with that of the DE, what re- 
mains to be clarified in future studies is whether the general- 
ized eigenstate thermalization scenario is valid or not for rij. 

In Fig.[l0j we present a study equivalent to the one in Fig. [9] 
but for the momentum distributions of SFs and HCBs. The 
first feature that is apparent in Figs. [Tol a) and HUl d) is the con- 
trast between the shape of momentum distribution of SFs and 
HCBs. For small values of k, the former resembles a Fermi 
sea and the latter resembles a bosonic system, both systems at 
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FIG. 10: (Color online.) [(a),(d)] Momentum distribution functions 
(L = 50), [(b),(e)] Am stat " DE , and [(c),(f)] a s ™, for SFs (left column) 
and HCBs (right column), in quenches to the delocalized regime, 
\i - —1.5 — s- Xf = 1.0. In [(b), (e)], the main panels depict 



Am GMEDE , while the insets depict Ar 



The legends for the 



ensembles considered are the same in all panels and are reported in 
(d). The values of <5gme and Sme are the same for all L's and are equal 
to the ones in Fig. [6] 



finite but low temperature [52||. However, the high occupation 
ofm k and ra\ in the tails makes it evident that the systems are 
not in thermal equilibrium, as can be seen by comparing them 
to the predictions of the ME. By contrast, the GME predic- 
tions for the momentum distributions closely agree with the 
DE results. This is better seen in Figs.fTOtb) andfTUle), where 
the differences between the GME and the DE predictions are 
seen to be small and decrease with increasing system size. The 
insets to the same panels show that Am MEDE is several times 
larger than Am GMEDE for the system sizes studied, and thai 
the fermionic Am ME_DE has a tendency to saturate to a finite 
value as L — >• oo, though the behavior of Am ME ~ DE for HCBs 
is less clear. 

The results for the scaling of a^f with increasing systerr 
size [Fig. [TUtc) and[T0Tf)l make apparent a fundamental dif- 
ference between the behavior of eigenstate expectation values 
of rh k and m k . For all ensembles, er^ at for SFs [Fig. [TUlc) 
is seen to be larger than the corresponding variance cr^f foi 
HCBs rFig. llOf f)!. and the former appears to saturate to a fi- 
nite value (more clearly seen within the ME than within the 
DE and the GME), while the latter appears to vanish (more 
clearly seen within the GME and the ME than within the DE) 
in the thermodynamic limit. This makes evident that local- 
ization of the single-particle fermionic eigenstates in momen- 
tum space has a clear consequence on the behavior of m{ 
in the many-body eigenstates, for which a^ 1 may be finite 
in the thermodynamic limit. By contrast, such a localization 
phenomenon in the single-particle basis appears not to have 
any effect on the many-body eigenstate expectation values of 



the momentum distribution of the HCBs, for which general- 
ized eigenstate thermalization can take place as ct^" appears 



{m b k ) DE 



(ml) gme 123H 



to vanish for all ensembles, and 

Again, we stress that even if the bosonic er^ E does vanish in 
the thermodynamic limit, we can infer that eigenstate thermal- 
ization does not occur for ml from the failure of thermaliza- 



tion to occur for this observable |33 



Quenches to the localized regime 

Density profiles obtained within the DE, GME, and ME af- 
ter a quench to the localized regime are shown in Fig. QTja), 
and the scalings of Art GME " DE ( An ME " DE ) with increasing sys- 
tem size are reported in the main panel (inset) in Fig. HTl b). 
Despite the fact that the site occupations fluctuate from site 
to site much more in Fig. [TTT a) than in Fig. Ufa), the GME 
results still closely agree with those of the DE, while the ME 
results do not. The scaling of An GME DE with increasing sys- 
tem size suggests that the differences between the predictions 
of the GME and the DE will vanish in the thermodynamic 
limit. The results for An ME ~ DE are less conclusive, although 
the values of this quantity obtained for the largest system sizes 
suggest a possibly tendency toward saturation. 

A clear difference between the behavior of the site occu- 
pations in the localized and delocalized regimes is seen in 
the fact that the variance of the eigenstate expectation values 
of this observable saturates in all three ensembles to a finite 
value in the former [Fig. [TTT c)l. whereas our results suggest 
that it vanishes in the latter [Fig.[9jc)], as the system size is 
increased. The saturation observed in Fig. [UJc) is a clear 
consequence of localization of the single-particle eigenstates 
in real space for A > A c . Finite values of the fermionic cr^ E 
in the delocalized regime and of a^ E in the localized one are 
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FIG. 1 1 : (Color online.) (a) Density profiles (L = 50), (b) An slatDE , 
and (c) a^ at in quenches to the localized regime, A/ = 1.0 — > Xf = 
4.0. In (b), the main panel depicts An GMME , while the inset depicts 
An ME DE . The legends for the ensembles considered are the same in 
all panels and are reported in (a). The values of <5gme and Sme are the 
same for all L's and are equal to the ones in Fig. [6] 
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FIG. 12: (Color online.) [(a),(d)] Momentum distribution functions 



(L = 50), [(b),(e)] Am sl 



, and [(c),(f)] aZ\ for SFs (left column) 



and HCBs (right column), in quenches to the localized regime, Xi = 
1.0 -> X F = 4.0. In [(b),(e)], the main panels depict Am GME " DE , 



while the insets depict Am 



The legends for the ensembles 



considered are the same in all panels and are reported in (d). The 
values of (5gme and Sme are the same for all L's and are equal to the 
ones in Fig. [6] 



physically relevant examples of the failure of the variance of 
few-body observables in the many-body eigenstates that con- 
stitute the ME to vanish in the thermodynamic limit, contrary 
to what is generally expected to occur for few-body observ- 
ables n. 

In Figs. fT2T a) and fT2l d). we show the momentum distribu- 
tion functions of SFs and HCBs in the DE, the GME, and the 
ME after a quench to the localized regime. In each ensemble 
the results for SFs and HCBs are barely distinguishable frorr 
each other, which we might intuitively attribute to localiza- 
tion undermining the particle statistics. As in the quenches tc 
the delocalized regime, the GME results closely follow those 
of the DE, and the main panels in Figs. [T2Tb) and [TSl e) in- 
dicate that Am GME DE decreases with increasing system size. 
Whereas we expect Am GME_DE to vanish for SFs in the ther- 
modynamic limit (as the DE and GGE predictions coincide 
for all one-body fermionic observables), for HCBs we expeci 
it to converge to a small but finite value, because of the failure 
of the GGE to describe the time averages of the HCB momen- 
tum distribution function observed in Ref. 113311 - In Figs.[T2la' 
and[T2Td). the ME results for the momentum distributions are 
clearly distinct from those of the DE, and that difference is ex- 
pected to remain in the thermodynamic limit both for SFs anc 
HCBs, as suggested by the scaling of Am ME ~ DE in the insets 
inFigs.[T2l'b)and[T2l;e). 

Results for the scaling of er^" with increasing system size 
are presented in Figs. [T2l c) and [T2l f). for SFs and HCBs, re- 
spectively. Once again, the results for the two particle species 
are very similar to each other. They are particularly inconclu- 



sive for c^ E , which is seen to increase with increasing system 
size for all systems except the two largest ones, for which it is 
seen to decrease. On the other hand, er° E and <r™ E decrease 
for all system sizes, though the results for the largest two sys- 
tem sizes suggest that they may saturate. This leaves open the 
question of whether af r f vanishes in the thermodynamic limit 
or whether it remains finite. What is clear from the fact that 
both the GE and the GGE fail to describe ml after relaxation 



M33I1 is that neither eigenstate thermalization nor generalized 
eigenstate thermalization take place in the HCB system in this 
regime. 



Quenches to the critical point 

For completeness, we present here results for the ensemble 
expectation values of observables in quenches to the critical 
point. We note that from the analyses of Sec.HlTland Ref. 0311 . 
we expect the critical regime to be most sensitive to finite size 
effects. 

In Fig. [13] we present results for the density profiles. They 
are in between those for quenches to the delocalized and lo- 
calized phases, as expected due to the finite and small nature 
of the lattice sizes studied. In particular, the average vari- 
ance er^f ' decreases with increasing system size but has a ten- 
dency to saturate, so much larger system sizes will be needed 
to resolve whether it vanishes in the thermodynamic limit 
(as it may in quenches to the delocalized phase) or whether 
it remains finite (as expected from the observed behavior in 
quenches to the localized phase). 

Figure [14] shows results for the momentum distribution 
functions of SFs (left columns) and HCBs (right columns). 
They are also intermediate between those obtained in the de- 
localized and localized regimes. The average variance er^" is 
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FIG. 13: (Color online.) (a) Density profiles (L = 50), (b) An slatDE , 
and (c) cr^ at in quenches to the critical point, Aj = —0.5 — > Xf = 



1.0. In (b), the main panel depicts An 



while the inset depicts 



An" 



The legends for the ensembles considered are the same in 



all panels and are reported in (a). The values of <5gme and <5me are the 
same for all L's and are equal to the ones in Fig. [6] 
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FIG. 14: (Color online.) [(a),(d)] Momentum distribution functions 
(L = 50), [(b),(e)] Am stat - DE , and [(c),(f)] cC', for SFs (left column) 
and HCBs (right column), in quenches to the critical point, Xi — 
-0.5 -»■ Af = 1.0. In [(b),(e)], the main panels depict Am GME " DE , 



while the insets depict Am 



The legends for the ensembles 



considered are the same in all panels and are reported in (d). The 
values of Sqme and Sme are the same for all L's and are equal to the 
ones in Fig. [6] 



larger for SFs than for HCBs, as in the delocalized phase. It 
can also be seen to decrease with increasing system size (as it 
does in the localized phase), but has a tendency to saturate (as 
it does in the delocalized phase), so much larger system sizes 
will be needed to resolve whether it vanishes or remains finite 
in the thermodynamic limit. 



V. SUMMARY 

We have studied the dynamics of the momentum distribu- 
tion function of noninteracting spinless fermions following 
quenches to the delocalized, localized, and critical regimes 
of a quasi-disordered lattice system. We found that although 
the time-averaged value of this observable agrees exactly with 
the predictions of the GGE in all three regimes, it does not 
exhibit relaxation after a quench in the delocalized phase. 
This is complementary to the failure of the on-site density of 
hard-core bosons (and therefore the noninteracting fermion 
system considered here) to equilibrate that was previously 
observed in Ref. [13311 - These behaviors can be understood 
in terms of localization of the single-particle eigenstates of 
the fermion model, in momentum space for the delocalized 
regime, and in real space for the localized regime, as discussed 
in Refs. I30ll34ll35ll . This analysis also helps us to understand 
the failure of the GGE to describe the momentum distribution 
functions of HCBs after relaxation in the localized regime ob- 



served in Ref. 0311 as a consequence of non-vanishing time 
fluctuations of one-particle fermionic correlations that persist 
in the thermodynamic limit. 

We found that in the delocalized and localized regimes, the 
SF observables that do exhibit relaxation to the GGE — the 
density in the former case, and the momentum distribution 
function in the latter — do so in a manner consistent with 
the Gaussian equilibration picture of Campos Venuti and Za- 
nardi 114011 . In quenches to the critical point of the Aubry- 
Andre model, we observed that both the density and the mo- 
mentum distribution function of SFs exhibit equilibration to 
the GGE, but that the decay of the time fluctuations of these 
quantities with increasing system size is slower than that pre- 
dicted by the conjecture of Ref. 14011 . 

We also studied the expectation values of one-body observ- 
ables in the many -body eigenstates of the SF and HCB Hamil- 
tonians, comparing results for the diagonal, microcanonical, 
and generalized microcanonical ensembles. Our observations 
indicate that the single-particle localization — in momentum 
space in the delocalized regime and in real space in the local- 
ized regime — that precludes relaxation of the corresponding 
observables to the GGE also leads to finite thermodynamic - 
limit variances of momentum and site occupations, respec- 
tively, in the many-body eigenstates of the SF Hamiltonian. 
By contrast, the failure of the GGE to describe the momentum 
distribution of HCBs after relaxation in the localized phase 
was not found to be associated with any corresponding satura- 
tion of the variance of momentum occupations in eigenstates 
of the HCB Hamiltonian with increasing system size. 

Because of finite-size effects, we did not find clear indica- 
tions of the behavior of the maximum difference between the 
expectation values of the density and momentum distribution 
in all the eigenstates contributing to the ensemble averages as 
the system sizes is increased. It would be particularly impor- 
tant to understand whether the values of observables in the 
individual eigenstates constituting the GME approach their 
average values in this ensemble with increasing system size, 
in order to clarify whether generalized eigenstate thermaliza- 
tion occurs in the delocalized regime — which would explain 
why the GGE works there for describing the momentum dis- 
tribution functions of HCBs after relaxation — and whether 
generalized eigenstate thermalization fails (as we expect) in 
the localized regime, where the GGE fails. It would also be 
interesting to see how the addition of nearest neighbor interac- 
tions 115711 . which break integrability, modify our findings for 
both SFs and HCBs. 
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